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Dynamical decoupling has been actively investigated since Viola first suggested using a pulse 
sequence to protect a qubit from decoherence. Since then, many schemes of dynamical decoupling 
have been proposed to achieve high-order suppression, both analytically and numerically. However, 
hitherto, there has not been a systematic framework to understand all existing uniform 7r-pulse 
dynamical decoupling schemes. In this report, we use the projection pulse sequences as basic 
building blocks and concatenation as a way to combine them. We derived a new concatenated- 
projection dynamical decoupling (CPDD), a framework in which we can systematically construct 
pulse sequences to achieve arbitrary high suppression order. All previously known uniform dynamical 
decoupling sequences using tt pulse can be fit into this framework. Understanding uniform dynamical 
decoupling as successive projections on the Hamiltonian will also give insights on how to invent new 
ways to construct better pulse sequences. 
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I. INTRODUCTION 

One of the major difficulties in the realization of quan¬ 
tum computing and quantum information processing is 
protecting the quantum state from decoherence. Quan¬ 
tum error correction protocols were developed to meet 
this challenge n For a review and recent progress of 
this area, see references B0. 

While quantum error correction can be regarded as 
a form of closed-loop control, dynamical decoupling has 
been proposed as a way to counteract the interaction 
between a quantum system and the environment by an 
open-loop control field. The idea of using pulse se¬ 
quences, to protect nuclear spins from classical decoher¬ 
ence, dates back to 1950, when the spin-echo method was 
found B . Since then, many pulse methods have been de¬ 
veloped in nuclear magnetic resonance spectroscopy Q. 
In 1998, it was first pointed out that a similar technique, 
periodical dynamical decoupling (FDD), can be applied 
to open quantum systems [9|. By using a control field, 
with duration shorter than the time scale of the environ¬ 
ment, dynamical decoupling can suppress the interaction 
between a qubit and the bath Q. Important theoret¬ 
ical progress was made by understanding the effect of 
dynamical decoupling as the result of a symmetrizing 
procedure over the group composed of all independent 
TT pulses [Eini. However, the finite switching time in 
real experimental conditions makes the symmetrizing im¬ 
perfect. Therefore, concatenated dynamical decoupling 
(CDD) was proposed to eliminate higher-order errors in 
the interaction Hamiltonian dill. 

In realistic experiments, the imperfection of pulses 
such as the rotation angle error and the finite width, 
must be taken into account. CDD is preferred over FDD 
in this case, due to the fact that concatenation not only 
suppresses the interaction with the environment but also 
the pulse errors to higher order d) d ■ 


Meanwhile, by optimizing the pulse interval to sup¬ 
press the low-frequency region of the noise spectrum, 
Uhrig dynamical decoupling (UDD) can achieve the same 
suppres sion order with fewer pulses, compared with CDD 
d) d ■ Although it has superior performance, UDD is 
very sensitive to pulse errors, due to the fact that it only 
uses single-axis rotations. However, to protect unknown 
states, uniform DD with multi-axis rotations can com¬ 
pensate errors due to its symmetric structure d> d-dl ■ 

Another development in dynamical decoupling is to use 
random pulses, instead of deterministic schemes, for suf¬ 
ficiently long sequences di d • Instead of trying to sup¬ 
press the interaction to arbitrarily high orders, random 
dynamical decoupling schemes improves the time depen¬ 
dence of the error accumulation from quadratic to linear 

dd- 

In this paper we only consider deterministic uniform 
dynamical decoupling for two reasons: (1) it is easy to im¬ 
plement in experiment, and (2) it is robust against pulse 
error compared to UDD. Besides analytical calculations, 
recently many other DD sequences have been found using 
genetic algorithms to optimize the suppression order [l^ , 
some of which even achieve same suppression order with 
fewer pulses than CDD. However, a unified understand¬ 
ing of all known DD schemes has not been developed 
until now. 

In this work, we propose concatenated-projection dy¬ 
namical decoupling (CFDD) to unify all known uniform 
DD schemes. Our framework gives a way to construct 
new pulse sequences and to calculate their suppression 
order. In Sec. im we describe the mathematical settings of 
the dynamical decoupling technique. In Sec. IHIl we first 
define our projection pulse sequence and explain its ef¬ 
fect as an ‘atomic’ projection. Then we use concatenated 
projections along different directions to construct more 
complex pulse sequences with arbitrary suppression order 
in Sec. m These results comprise the two cornerstones 
of our theory of CFDD. In Sec.|V]we formally introduce 
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CPDD by defining the CPDD equivalence classes, which 
are specified by three integers. In subsection A, a series of 
properties of CPDD are developed; in subsection B we de¬ 
sign a deterministic scheme to construct optimized pulse 
sequence for each suppression order. A table of known 
DD schemes is given as well to show how these known 
schemes fit into our CPDD framework. In Sec. ED we 
discuss why, intuitively, some CPDD sequences are su¬ 
perior than CDD. We also point out a typo in Ref, [l^ . 
which is easily detected within the framework of CPDD. 
We present our conclusion in Sec. I VIII 


The first two terms are 

Hio) = 1 ° H{t)dt, 

A Jo 

rdh [ ^ dh[H(h),H{h)], 

Jo Jo 

We make the separation, 

H = Ti-s ^ Hb Hen, 


( 6 ) 

( 7 ) 


II. DYNAMICAL DECOUPLING SETTINGS 

We consider a qubit system S coupled to an arbitrary 
bath B, which forms a closed system on the Hilbert 
space Hs ®Hb- The overall Hamiltonian can be written 
in the form, 

Ho = Hs 0 IIB -I- Is ® Hb + Hsb, (1) 

where 1 is the identity operator, Hs is the pure system 
Hamiltonian, Hb is the pure environment Hamiltonian, 
and Hsb is the interaction between the qubit and the 
bath. In the following, we will assume the interaction 
takes the general linear form. 


in which Herr is the part responsible for undesired evolu¬ 
tion of the qubit. 

In the context of dynamical decoupling, the only mea¬ 
sure of the performance of a pulse sequence is the sup¬ 
pression order. The suppression order N is defined 
through the following equation, 

OO 

Herr= ^ ( 8 ) 

which means the corresponding pulse sequence achieves 
Nth Qj.(jgj. decoupling. 

III. PROJECTION PULSE SEQUENCE 


Hsb = cTx'S’+ ay igi By + a,, B,,, ( 2 ) 

where B ^, By and B^ are abitary operator on Hs. 

Decoherence can be suppressed by adding a control 
field solely to the system, i.e., Hc(t) (8i 1 _b- If the control 
field is a series of pulses, then the suppression effect can 
be understood as a symmetrizing procedure M- In the 
toggling frame, the Hamiltonian is transformed into. 


Now consider that we use K uniform ideal tt pulses 
as the control field. Thus the field takes the form 
Hc{t) = - tj), where tj = jra, Td is the 

pulse interval and aj G {1, Uy, Uz}. In the following, 
we use Pk---Pi, Pk G {/, A, V, Z} to represent corre¬ 
sponding pulse sequence. In the limit of Td 0, we have 
a continuous pulse sequence, and only the zeroth-order 
term in the Magnus expansion survives. 


H{t) = Ul{t)HoUe{t), (3) 

where Uc{t) = Texp{—i J* Hc{t)dt} is the evolution 
operator of the control field. Here T is time ordering 
operator. The evolution of the state vector is governed 
by the evolution operator, 


h(o) = - [ ^ UUt)HoUeit)dt. (9) 

Tc Jo 

Also, since the pulses are ideal, Uc{t) is a piece-wise 
constant function. The first-order average Hamiltonian 
reduces to 


C/(t) = Texp I — i J H{t)dt^. 


( 4 ) 


( 10 ) 

i=i 


If the sequence length is Tc, we can define an average 
Hamiltonian H{tc) by U(tc) = exp^—iHrc). By expand¬ 
ing U{tc), we can collect terms according to different 
order of Tc, 

Texp { - * ^ (5) 


If we chose pulses such that Uc{tj) go through each ele¬ 
ment of certain group such as G = {I, X, Y, Z}, Eq. (fTOll 
is just a symmetrizing procedure which projects Ho onto 
the commutant of the group algebra [1^ . 

A more intuitive way to view the effect of pulse se¬ 
quence is to look at them as combination of basic projec¬ 
tions [l^. The simplest pulse sequence is Pi Pi, which is 
similar to the CPMG pulse sequence BUI. Hereafter we 
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call it the projection pulse sequence and use the notation 
Pi = PiPi- To explain its projecting effect, we consider 
the spin-boson model, which induces the longitudinal de¬ 
cay of the spin. The interaction takes the form, 

Hsb = <S'bk +gl(T~ (Sibl), ( 11 ) 

k 

where bk is the annihilation operator of a photon with 
momentum k, and gk is the coupling strength between 
photon with mode k and the spin. Here is the creation 
(annihilation) operator, = Ux Pioy. 

If we apply pz = ZZ to the system, then Uc{tj) G 
Q = {1, Z}. Using Eq. (|T0)) . the interaction term will be 
completely removed in the continuous limit due to the 
fact CFz<Jx(y)0'z = Therefore, geometrically the 

pulse sequence Pz projects the Hamiltonian along the z 
direction. 

However, in real experimental conditions, there is an 
upper limit of the pulse switching rate; thus is finite. 
Although the Magnus expansion is still valid, so long as 
TcLOc << 1, where Wc is the cut-off frequency of the bath, 
the projection is not exact anymore due to higher-order 
corrections. 

Theorem 1 .■ Assume the interaction between a single 
qubit and bath takes the form of Eq. (0). After applying 
the projection pulse sequence pj (j = x,y, z) with pulse 
interval Td, the zeroth order error Hamiltonian defined 
in Eq. & and Eq. m is given by, 

Hi°}=7TfHo = a,<^B,. ( 12 ) 

Thus the full error average Hamiltonian is given by, 

-Herr = CTj (g) [Bj + + X! ® 

iJ-j 

(13) 

where we use to represent the mapping from Hq to 

Herr that is induccd by the projection pulse sequence pj. 
The symbol T represents the directions orthogonal to di¬ 
rection j and B^^^ ~ k^'^oljd). Here is some combi¬ 
nation of commutators of the bath operators with dimen¬ 
sion of [H]^. 

This projection point of view gives a geometrical and in¬ 
tuitive way to understand the effect of the pulse sequence 
PjPj. The full proof is given in Appendix A. 

In summary, we have shown that the effect of the pro¬ 
jection pulse sequence pi = Pi Pi is to project the Hamil¬ 
tonian along i direction up to first order. 

IV. CONCATENATION OF CYCLIC PULSE 
SEQUENCES AS SUCCESSIVE PROJECTIONS 

The higher-order terms remaining in the -Herr, after ap¬ 
plying the projections, will coherently add up with time. 


To achieve higher-order suppression, we need to project 
the Hamiltonian along different directions successfully. 
We will show in this section that the effect of the concate¬ 
nation of cyclic pulse sequences is to apply successively 
the projections induced by each pulse sequence. 

A pulse sequence is called cyclic when the generated 
evolution operator is periodic with period Tc up to a 
phase factor, 

Ucinr,) = (14) 

where n = 0,1, 2 ... and (j) is an arbitrary phase. 

An equivalent definition of a cyclic pulse sequence is 
that the product of all the pulses is equal to one, up to 
an arbitrary phase, 

K 

l[n = PH, (15) 

which follows directly from Eq. (II4p when n = 0. We will 
see that this property is necessary for the proof of the 
equivalence between concatenation and successive pro¬ 
jections. 

Another useful property of the cyclic pulse sequence 
is that the concatenation of two cyclic pulse sequences is 
still cyclic. The proof is straightforward and we include it 
in the Appendix A. Having the definition of cyclic pulse 
sequence, we now prove the second basic theorem of our 
CPDD scheme. 

Lemma 1 .• Consider two pulse sequences A and B, 
Pf^.'.Pl, where i = A, B, with the same pulse interval. 
The first pulse sequence A = PjPj (j G {x,y,z}), which 
is a projection pulse sequence, and sequence B is concate¬ 
nated from multiple projection pulse sequences. A third 
pulse sequence C is constructed by concatenating A and 
B, C = A[B] = PjBPjB. The following relationship is 
true, 

TTc = ( 16 ) 

where the mapping 7r° induced by applying the pulse se¬ 
quence i is defined in TheoremUi 

Lemma [1] is the theoretical cornerstone of this work. It 
explains why concatenation can increase the suppression 
order, which is not so obvious. The cyclic properties and 
the changeability of different pulses (up to an irrelevant 
phase factor) are necessary for the proof. The full proof 
is included in the Appendix A. 

V. CONCATENATED PROJECTIONS 
DYNAMICAL DECOUPLING 

What really distinguishes our work from the CDD 
scheme is that we chose the projection pulse sequence as 
the basic element of concatenation. Motivated by Theo¬ 
rems [T] and [5J we define our concatenated-projection dy¬ 
namical decoupling (CPDD) as a new way to construct 
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pulse sequences by applying projections along different 
directions successively. Since each projection kills the in¬ 
teraction terms orthogonal to it by one more order, by 
appropriately combining different projections, our CPDD 
can achieve arbitrarily high suppression order. 

Definition 1 .• A CPDD pulse sequence is specified by 
an ordered series ipf, In-Ij H ■ It constructed by con¬ 
catenating N projection pulse sequences successively, A = 
PiN[ViN-A‘’\Pii ]•••] ^ where ij € {x,y,z} andl <j<N. 

N 

The suppressing effect of the CPDD sequence on the 
Hamiltonian follows immediately from the combination 
of the effects of projections and concatenation. 

Theorem 2 ; Consider a CPDD pulse sequence A spec¬ 
ified by i]\j,ipf-i, ...,ii. After applying pulse sequence A, 
the average interaction Hamiltonian is given by 

^ (17) 

i=x,y,z 

where 

di = 'y ^ Uj , (18) 

j_Li 

and Uj is the number of pj sequences. 

Proof : Repeatedly using Lemma [U the leading order 
of the error average Hamiltonian after applying sequence 
A is given by 


by each of them are He?r = J2i=x y ® and 

Sirl = y z^i ® define A and A' to be 

equivalent to each other 

A - A', (22) 

if the leading order of the average error Hamiltonians are 
the same along each direction, namely ni = 

It can be easily proved that the relationship defined above 
satisfies the three properties of an equivalence relation¬ 
ship. 

Therefore, for a CPDD sequence specified by sequence 
a = jat, *Ar-i,..., * 1 , all CPDD sequences specified by 
a’s permutations A{iiq,ij^-i, form an equivalence 

class. By the virtue of the equivalence class, only three 
numbers nx,ny,nz are needed to completely specify a 
CPDD class. 

Definition 3 .• A CPDD class is defined as an equiva¬ 
lence class with equivalence relationship defined in Def- 
inition\^ specified by three integers, {nx,ny,nz}. The 
structure of the pulse sequence can be generated by con¬ 
catenating all Hi Pi sequences (i = x,y,z) in arbitrary 
order. 

From the definition of CPDD and Theorem [5J we derive 
a series of properties satisfied by CPDD sequences and 
their equivalence classes. 

A. Properties of CPDD 


( 0 ) „ ( 0 ) ( 0 ) ( 0 ) jx 


(19) 


From Theorem [U each projection remove the first order 
term in the perpendicular direction. Therefore, 

Hi?r= E (20) 

i—x,y,z 


where 

di = 'y2nj- ( 21 ) 

j±i 


Due to the way concatenation connects two pulse se¬ 
quences, we can derive two properties that the structure 
of each CPDD sequence must satisfy. 

1. For an arbitrary CPDD pulse sequence, each 
odd site has the same kind of tt pulse. 

We prove this by induction. Consider a pulse sequence 
A = PkPk-i---Pi, which is concatenated from N pro¬ 
jection pulse sequences. Pulse sequences are concate¬ 
nated from the first n (1 < n < N) of them. The follow¬ 
ing relations are satisfied, 

An=Pi[An-l], (23) 


From Theorem [5] we can see that the effect of ni 
PiS is to suppress the error Hamiltonian along the i 
direction to nAh order. From Eq. (ED we also notice 
that the order of how different projection pulse sequences 
is concatenated does not affect the leading order of the 
average Hamiltonian along each direction. Therefore we 
can define an equivalence relationship between different 
CPDD pulse sequences. 


where pi is the nth projection pulse sequence. Assume 
for subsequence the pulses are the same for each 

odd site, 

An-l = P2n-2pb.p2Po- (24) 

Let’s examine the pulse sequence 
An — Pi [kin—l] 

= (P,P2„-2)Po-P2Po(P.P2n-2)Po-P2Po. (25) 


Definition 2 .• Consider two pulse sequences A and Therefore, the fact that all the pulses at each odd site 
A'. The leading order of the error Hamiltonians induced are the same still holds. 
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Since A 2 = Pi[pj] = {PiPj)Pj{PiPj)Pj, which also has 
the same kind of pulse on its odd sites, by induction we 
have proved that for each odd sites of Ajv the pulses are 
the same: 

^ 2 m+i=^b, m = 0,l,... (26) 


2. For an arbitrary CPDD pulse sequence, the 
first half and the second half subsequences are the same. 
Again this property follows from the definition of 
concatenation. Using Eq. (l23ll ioi n = N we have 

A = p,„[A„_i]. (27) 

Assume A„_i = P 2 n- 2 P 2 N- 3 ---Pi, we have 

^ = {Pi„P2N-2)P2N-3---Pl{Pi-nP2N-2)P2N-3---Pl- 

(28) 

Obviously sequence A is composed the same two copies 
of Aat-i, or more precisely 

Pm = Pm-N/2, m = l,2,...,N. (29) 


3. For CPDD class {u^,, rij,, n^}, the number of 
pulses or sequence length K is given by 


The proof is straightforward. At first, each basic 
projection is induced by the two same pulses. Secondly, 
the length of two concatenated pulse sequences, A 
and B, is equal to the product of the length of each 
two, Ka[b] = KaKb- Therefore, for a pulse sequence 
composed of rii pairs of {Pi, Pi), the total pulse number 
K is given by Eq. (1501) . 

4. For the CPDD class {nx,ny,nz}, the suppres¬ 
sion order achieved is given by, 

= min [Uy + Uz, Ux + n^, Uj, + Uy}. (31) 

From Theorem [21 the leading order of the error Hamil¬ 
tonian induced by any pulse sequence in the CPDD class 
{na;,ny,n„} is given by 

E (32) 

i—x,y,z 

where di = Since the suppression order N 

is defined as the leading order of H, N = min {dx, dy, dz}. 

From the expression of suppression order, Eq. ®, 
we can see that simply increasing pulse numbers (num¬ 
ber of projections) does not necessarily increase the 


suppression order. Actually, in the framework of CPDD, 
only for pulse sequences with certain pulse numbers can 
the suppression order increase. 

5. For a given suppression order N, the minimum 
number of pulses, A'min, required to achieve such sup¬ 
pression order is given by 


l0g2{KnAn) = i[3iV + ^(1 - 1®"^)], (33) 

where 1®'^ = 1 © 1 © ... © 1. 

^ ^ > 

N 

To achieve suppression order N, SN terms in the 
interaction Hamiltonian need to be eliminated due to 
the form of Hsb, Eq. ©■ However, each basic pro¬ 
jection TTi requires two pulses, which implies that only 
an even number of terms can be eliminated for a given 
CPDD sequence. Therefore, we need to add one more 
pulse depending on whether N is odd or not. Divid¬ 
ing the total number of eliminated terms by two we have, 


^ „, = i[3iV+1(1 +(-1)^+1]. (34) 

i—x,y,z 

Using the results of property 3, we have Eq. (1551) . The 
mysterious series 4,8, 32,64, 256... was first found by a 
genetic algorithm in Ref. [l9l| which is now understood 
within the framework of our unifying CPDD. 


B. Optimized uniform dynamical decoupling 

The pulse sequences corresponding to the pulse num¬ 
ber in Eq. (|33|) uses the minimum number of pulses at 
each suppression order. This optimized uniform dynami¬ 
cal decoupling (OUDD) scheme can be represented using 
the CPDD indexes as 

OUDDfc : l{fc-l®^fc + l®'=,fc + l®'=}. (35) 

The suppression order of OUDD^ is Nk = k and the 
sequence length is given by ATmin in Eq. (I55|) . 

As we can see, some of the pulse sequences are 
particular levels of CDD; and GA8;. To compare with 
OUDD and other known DD schemes, we list the 
corresponding CPDD indexes of known DD schemes in 
Table |T] below. 
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TABLE I: Known DD schemes represented as 
CPDD 


\tIx , Tiy , Tlz } 

Name 

Pulse sequence 

K 

N 

{0,0,1} 

Projection 

P^P^ 

2 

0 

{0,1,1} 

PDD(CDDi) 

PiPjPiPj 

4 

1 

{1,1,1} 

GA8a 

I Pi Pj Pi I Pi Pj Pi 

8 

2 


CDD; 

CDD[CDD/i] 

4 ' 

1 

{1,1,1} 

GA8, 

GA8a[GA8,i] 

8' 

21 


VI. DISCUSSION 

Although CDD also relies on concatenation, the 
fact that our CPDD uses basic projections as building 
blocks makes finding more efficient pulse sequences 
possible. To make this clear, we consider the CDD from 
the view point of our CPDD. CDDi can be consid¬ 
ered as the concatenation of two different projections, 
CDDi = X(YY)X(YY) = ZYZY [II]. Therefore the 
effect of CDDi is to apply projections along the y and x 
direction successively, 


= [(^x (Z B^ +ay Z) +az® B^^] 

= ax® B^^'> +ay ® B W + a^ ® B^^'> 

- fc(i)o(T]). (36) 

As we can see, CDDi completely removes the zeroth or¬ 
der interaction terms, thus achieving suppression order 

.^CDDi = 1- 

Now consider CDD 2 , which is the concatenation of 
two CDDi sequences: we write explicitly the process of 
the successive projections , 


t^CDD^Hsb = 7r®7r(°V(°)7r),°)iLo 

= [ax ® + ay® B^^ +a^® Bf">\ 

= ax® B^^'> +ay® B^^'> +az® B^^^ 

- (37) 

As we can see from above, a total of eight eliminations 
(each projection pulse sequence eliminates two terms 
in the orthogonal directions) are used to completely 
remove the first two orders of the interaction Hsb- 
However, the two additional eliminations of does 
not contribute to further increasing of the suppression 
order. To avoid this, we consider projecting along each 
direction exactly once, namely TTxT^yT^z, which belong 
to the CPDD class {1,1,1}. Translating the projection 
back to corresponding pulse sequence according to the 
rule of concatenation, we have 


: Px[Py[Pz]] 

= Px[Y{ZZ)Y{ZZ)] 

= X{XZXZ)X{XZXZ) 

= IZXZIZXZ, (38) 

which only uses eight pulses and six projections. This 
was first found by a genetic algorithm and called GA8a 
in Ref. [T^ . 


To achieve suppression order of V = 2, CDD 2 needs 
16 pulses while GA8a sequences only requires 8 pulses. 
This efficiency of using pulses comes from the very fact 
that GA8a uses basic projections tt^ as building blocks 
while CDD 2 uses composed projections TTiitj {i ^ j) as 
building blocks. 

In Ref. [ 1 ^ , the author claimed to find another 8 
pulse sequence GXSh=Z{XYXY)Z{XYXY) which also 
achieved suppression order oi N = 2. From the structure 
of GA8b, we know the projections induced by it is, 


TT 


(0) 

G.486 


^ ( 0 ) ( 0 ) ( 0 ) 
z "■ z "y ; 


(39) 


which belong to CPDD class {0,1, 2}. Using the results 
from Eq. m, the suppression order of GA8b is equal to 
1. Therefore the claim in Ref. [l^ is typo. 


To double check our results, we also use the multi- 
precision package mpmath [28l| to compute the suppres¬ 
sion order of both GA8a and GA8b for a 5-spin model 
with random coupling constants. Here the distance D 
is defined as the distance between an actual evolution 
operator and the unit operator (T^ . 

D(U, Is) = ^i-^lirilT., (40) 

where d-u is the dimension of the Hilbert space, P = 
Trs[17] and || • Htj- represents the trace norm. An upper 
bound of D can be calculated (T^. 


D < 0[t^+\ (41) 


Therefore, we can extract the suppression order by 
plotting D versus in the log-log diagram. As we can 
see in Figure [U GAgo achieves higher suppression order 
than GAgb which is in agreement with the argument 
from the view point of our CPDD. 


VII. CONCLUSION 

We have developed tour CPDD schemes in which 
CPDD sequences are concatenated from different 
projection pulse sequences. We also define CPDD 
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FIG. 1: Comparing the suppression order of GA8a 
and GASb. D is the distance between actual evolution 
operator and the unit operator, defined in 0 , and 
Td is the pulse interval. The suppression order N is 
defined by the relation D ~ We consider 

the parameters J in the range of Jrd € 10“^], 

where J is the norm of the interaction Hamiltonian. 


equivalence classes as the set of pulse sequences that 
have the same leading order along each direction in the 
Magnus expansion of the average Hamiltonian. Based 
upon the definition of our CPDD pulse sequences, we 
prove a series of properties about the structure of the 
sequences that must hold for CPDD. We also give 
a formula to calculate suppression order for a given 
CPDD class specified by {nx,ny,nz}- We propose the 
optimal uniform DD sequence given in Eq. p5|) use in 
experiments, since each of the sequences in the series 
achieves its suppression order using minimum pulse 
numbers. Although some of OUDD are already known 
DD sequences, our CPDD framework gives a unifying 
and consistent way to both understand and construct all 
of them. 

The main advantage of using UDD is that the pulse 
number needed scales linearly with the suppression 
order, Nuod ^ Cl(iC), which is much more efficient 
than the exponentially dependence in CPDD. However, 
UDD is subject to several difficulties. Firstly it is 
valid o nly for environmental spectrum with a hard 
cut-off [ld . l^ , and secondly it is very sensitive to pulse 
errors |22 h^ . However, for our CPDD, especially 
the OUDD class, some pulse sequences have rotation 
symmetry thus making them robust against pulse 
error. Although we have not given a rigorous bound 
analysis for CPDD here, the results and the calculation 
should be similar to those in Ref. [0|: the distance 
between the actual state and the desired state goes to 
zero as the concatenation level goes to infinity. More 


importantly, pulse errors are suppressed along with se¬ 
ries of concatenations as long as the error is not too large. 

Moreover, although the suppression order of CPDD 
is still exponentially dependent on the pulse number, 
.Mcpdd 0 {\og 2 K), the reduction of the number of 
pulses for the same suppression order is exponentially 
large compared to the case of CDD, iVcDD ~ Oilog^K). 

Given the structure of our CPDD, one can see 
that two directions to try to find more efficient pulse 
sequences are worth exploring. Firstly, to find basic 
sequences, which could achieve higher N to K ratio than 
projection pulse sequences. The concatenated UDD 
Hi, [13 obviously belongs to this class. However, the 
applicability of Lemma 1 in the context of non-uniform 
DD is questionable and requires a new proof. Secondly, 
we seek to find new ways to combine two pulse sequences 
or two projections. Since the fact that exponentially 
large pulse numbers are needed to achieve high suppres¬ 
sion order resulting from the concatenation , new ideas 
in this direction may greatly improve the efficiency of 
uniform dynamical decoupling. 
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Appendix A 


In this appendix we include the proof of Theorem 
[T] and Lemma [1] which are the foundations of our 
CPDD scheme. A proof about a property of cyclic pulse 
sequences is also included. 

Proof of Theorem m Without loss of generality, 
we consider the pulse sequence is ZZ, where Z = —ia^. 
The transformed Hamiltonian H (t) is given by 


Hit) = 


Hi, 0 <t <Td 
H 2 , Td<t <2Td , 


(Al) 


where 

Hi = Bn <7x Bx + (Jy<Z By + az® Bz, (A2) 
H 2 = "I® Bn - (Tx® Bx - (Jy® By + az® Bz- (A3) 

After applying the pulse sequence PjPj, the first and sec¬ 
ond order expansion of the average Hamiltonian is given 







by Eq. dH), 


i^l 

(A4) 

HP^ = Y^eY.lP^K3]- 

(AS) 


i<j 


Using Eq. we have 

H^°'> = 1 (g) Bo + cr^ (g) Bz, (A6) 

= Td[Bo,Bx]^ax 

+Td[Bo, By] ® CTy + 2Td[By, Bx] (g) tJz. (A7) 

Since the commutator between bath operators is not zero 
in general, we have 

Herr = (^ z ® [^2 + k‘'H o{Td)\ + ^ CTi (g) [B^^^ + /cf ^o(tJ)], 

lA-z 

(AS) 

where B^^'^ = Td[Bid,Bi\. The same calculation gives 
similar results for j = a:,y. Q.E.D. 

Proof of Lemma [1] : The concatenated sequence C is 
given by 

{PrP^)Pi.^...P^{P^Pi)...P^, (A9) 

where the bracket means the there is no free evolution 

in between the two pulses inside the bracket. Since the 
projection pulse sequence is cyclic by definition, and se¬ 
quence B is also cyclic by Theorem O 

K 

(AlO) 

The 2Kb evolution operator of the control field, is 

given by, 

=n p? (All) 

j<m 


where we have used the commutativity of Pauli matrices 
and the cyclic property Eq. (lAlOl) of pulse sequence B. 

Now adding the action of and 
together, we have 

uZ^HsbUZ + 

= U^^HsBUi + Ug^P}HsBP^U^ 

= U^\2T:fHsB)U^ (A14) 

If m = Kb, we have and U^k^, which are 

= p^pk n pf 

3 <Kb 

= p. Xi p’’ 

j< = KB 

= “‘^PiVK, (A15) 

u?Ks = a( n pf)p{ n pf) 

3 <=Kb 3 <=Kb 

= (A16) 

Now using Eqs. I|A14[ lAlSi IA16p . the first order of aver¬ 
age Hamiltonian after applying sequence C is given by 

, 2 Kb 

^ E U^^PsbU^ 

m—1 

-j Kb 

m—1 

Kb ^ Ka 

= E 11™' ( 2 E Ui 

m—1 1—1 

= '^^b'^^aHsb (A17) 

Q.E.D. 

Theorem 3 The concatenation of two cyclic pulse se¬ 
quences is still cyclic. 


To construct and 7rg\ we group and 
together(m < Kb). 

For 1 < m < Kb, 

n™ = n pf=p^^ (A12) 

j<=m 

and, 

Kb 

= (n p,'’)«(n^^'’)(n pf) 

j<^—rn j'>'xn jd—m 

= p*p, n pf 

j<=m 

— c 2 J 


Proof of Theorem [3] : Consider two cyclic pulse se¬ 
quences A and B. From [Eq. (|T5|) ] we have 



II 

(T) 

(A18) 


Kb 

Y[Pf = 

2=1 

(A19) 

The pulse sequence C is constructed by concatenating A 
and B, thus 

C = A[B] 

■■PKB)P2^iPF-PKB)-PtiPl- 

-Pkb) 

(A20) 


(A13) 
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Therefore the product of all pulses of sequence C is 






Kb Kb 

i=l i=l 

pA i4>B pA i4>B 

Ka 

e^KB<PB JJ- pA 

(A21) 


where we define 4>c = Kb(I)b + 4>a- Therefore, the pulse 
sequence C=A[B\ is also cyclic. Q.E.D. 
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